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Abstract 

We present a barrier potential with bound states that is exactly solvable and 
determine the eigenfunctions and eigenvalues of the Hamiltonian. The equi- 
librium density matrix of a particle moving at temperature T in this nonlinear 
barrier potential field is determined. The exact density matrix is compared 
with the result of the path integral approach in the semiclassical approxima- 
tion. For opaque barriers the simple semiclassical approximation is found to 
be sufficient at high temperatures while at low temperatures the fluctuation 
paths may have a caustic depending on temperature and endpoints. Near 
the caustics the divergence of the simple semiclassical approximation of the 
density matrix is removed by a nonlinear fluctuation potential. For opaque 
barriers the improved semiclassical approximation is again in agreement with 
the exact result. In particular, bound states and the form of resonance states 

are described accurately by the semiclassical approach. 
PACS numbers: 82.20.Kh, 03.65.Sq, 05.30.-d 
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I. INTRODUCTION 



The semiclassical approximation to quantum mechanical problems is useful in a large 
variety of contexts. In chemistry and physics this approach becomes more prevalent for 
many reasons. For instance, semiclassical methods are efficient for the calculation of highly 
excited states, for which direct quantum mechanical calculations become difficult. The 
semiclassical approach also offers conceptual insights into the dynamics of many systems 
that are not easily extracted directly from the quantum mechanical treatment Finally, 
physical quantities of systems with barrier potentials can be calculated in the semiclassical 
limit, that corresponds to a large barrier height or a large barrier width 0. Here we 
study a system with a nonlinear potential field that has a barrier and an adjacent well. 
For the specific potential field considered, the classical equation of motion as well as the 
corresponding Schrodinger equation can be solved exactly. Therefore, the model can be 
used to test the quality of approximations to which one has to resort for most realistic 
potential fields. 

The imaginary time path integral approach provides a consistent scheme for the semi- 
classical approximation of several quantities in statistical mechanics P|. It has been noted 
already 25 years ago [^] that classical paths methods offer facile techniques to determine the 
semiclassical approximation of the equilibrium density matrix or the partition function. In 
the last few years the formulation of quantum statistical mechanics based on the centroid 
density, originally suggested by Feynman, was studied extensively p. Although this for- 
mulation allows for highly accurate approximations for the partition function, the centroid 
density has the same physical interpretation as the density matrix only in the high temper- 
ature limit. The simple semiclassical approximation for the density matrix studied in the 
past becomes exact at high temperature, however for coordinates near critical values di- 
vergencies arise when the temperature is lowered 0. Well-known as the problem of caustics, 
these divergencies refiect the fact that new classical paths become available. In the region, 
where these paths emerge, one has to improve upon the simple semiclassical approximation 
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and evaluate non-Gaussian fluctuation integrals 0. 

In this article we compare the exact equilibrium density matrix with the improved semi- 
classical approximation. First, in section 2 we calculate the exact equilibrium density matrix 
using analytic solutions of the Schrodinger equation. The simple semiclassical approximation 
is evaluated in section 3. In section 4 we present an improved semiclassical approximation 
for the density matrix which remains valid near the caustics. Section 5 presents semiclassical 
results for bound and resonance states. 



II. EXACT EQUILIBRIUM DENSITY MATRIX 

We consider the thermal equilibrium state of a quantum system with the asymmetric 
barrier potential 

, , A + 5sinh(X/Lo) 

cosh-(A7L:) ■ « 

where A and B are dimensionless coefficients, Lq denotes the typical range of the potential 
and the energy scale. The Schrodinger equation with potential (|1|) will be seen to be 
exactly solvable. In particular, for i? = we recover the symmetric Eckart potential 

The antisymmetric part of the potential proportional to sinh(X/Lo)/cosh^(X/Lo) is 
different from the tanh(X/Lo) term familiar from the work by Rosen and Morse This 
leads to important qualitative differences in the shapes of these exactly solvable potentials. 
As seen from Fig. 1, the potential (|l]) always approaches the same limiting value for X — *■ 
±oo, and for i? 7^ it displays a barrier with an adjacent well. The fact that the potential 
is a genuine scattering potential decaying exponentially for large |X| and the barrier- well 
form distinguish (||) from other exactly solvable smooth potentials. 



A. The Schrodinger Equation 



The Schrodinger equation of a quantum mechanical particle moving in the potential (|Tp 
may be written as 
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+ V{x) 



with the scaled potential 



V{x) 



a + 6sinh(x) 
2 cosh^fx) 



In (|^) and we have introduced the dimensionless quantities 

X = X/Lo e = mLlE/h^, 

a = 2mUoLlA/h'^, b = 2mU oLlB / . 

For later purpose we also introduce the dimensionless inverse temperature 



(2) 



(3) 



(4) 



(3 = h^/mLlk^T, 



(5) 



where k-Q is the Boltzmann constant |T0(| . 

With the transformation y = sinh(x) the Schrodinger equation takes the form 

2e a + by 



y 



l + y- 



^(y) = 0. 



(6) 



.1+1/2 (l + y2)2_ 

This second order differential equation has three singular points at i, —i, and oo. It can be 
reduced to its canonical form with the ansatz 



(7) 



where 



s = Vl - 4a + iAb = s' + is", 



(8) 



with the real and imaginary parts given by 

^ '[l-Aaf + iAb)' 
1 



1 - 4a 1 2 



s" = sign(6)^^[(l-4a)2 + (46)2 



1 - 4a 1 2 



(9) 



Inserting (|^ into (|^), we obtain with the substitution 
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a hypergeometric differential equation for F{z) 



z{l - z)F"{z) + [7 - (a + a* + l)z\F'{z) - aa*F{z) = 0. 



Here, the coefficients a and 7 are given by 



1 s' r- 



(11) 



(12) 



A solution of (|TT]) is given by the hypergeometric function 2-^"! («,«*; 7; 2^) |T2|- '^^^ 
corresponding solution of equation (^) reads 



^(2) = [i22]3+4^ [i2{z 2Fi{a,a*]r,z). 



(13) 



This solution is defined in the complex plane for \z\ < 1, i.e. for coordinates < arcosh(2). 
Beyond this region the solution is analytically continued by means of the formula 



w( * \ r(7)r(a*-a) „ 1, 

i [a*)L (7 — a) z 



r(7)r(a-a*) 



iFi{a*, 1 — 7 + a*; 1 — a + a* 



1, 



(14) 



r(a)r(7-a*)' ' 

This way we have found an exact solution of the Schrodinger equation defined in the whole 
domain of x . Of course, physically acceptable solutions have to satisfy proper boundary 
conditions. 

Let us first consider the case e > 0, i.e. a typical scattering situation. Since the potential 
(§) decays exponentially for large the wave function \E'(x) should reduce to plane wave 
states at large coordinates. For x ±00 the argument I/2; in ( p!4D vanishes and one can 
show that the asymptotic behavior of the solution (|13D is indeed of the form 



\E'(x ^ ±00) = A± exp 



-iV2e(±x) 



+ B± exp iV2e{±x) 



(15) 



with energy dependent amplitudes 



r(7)r(a* -a) f 1 + s' r, . ^ 



r(Q;*)r(7 - a) 



71 



± 



71 



±2i ln(2) 



and 



T{^)T{a-a*) i 1 + s' 



ln(2) T ^l 



71 



-±2iln(2) 



r(a)r(7-a*) [ 2 
Accordingly, the complex conjugate solution ^*{x) behaves asymptotically as 

±oo) = A*^ exp (±iV2€x^ + exp (^p^V^a;) , 

We superimpose these two solutions to wave functions 

^±{x) = u±'^{x) + v±'^*{x) 

with the asymptotic behavior 

^±{x — > ±oo) = exp (±iV2ex^ , 



describing outgoing scattering waves. The coefficients u± and v± are found to read 



v± = 



A, 



|5±|2- |A±|2 

The asymptotic behavior of ^±{x) for x — > ^oo is then given by 

^±{x — > ^oo) = — ^ exp (±iy/2ex) + ^^^J^^ exp (^iV2ex\ , 
t±{e) ^ ' i±(e) ^ ' 

where the transmission and reflection amplitudes t±{e) and r±(e) are given by 

t±(e) 



\B±?-\AA^ 



31B^ 



BIB^ - A*^A± 



The amplitudes t+ and r+ can be expressed in terms of r_ and t_ by 
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Clearly, t±{e)^±{x) describe the stationary scattering of particles incident from the far left 
and right, respectively, with energy e and unit amplitude. Inserting the explicit result for 
the amplitude factors A± and B±, the transmission probability T(e) = |t±(e)p is found to 
read 



T(e) = exp (vrs") 



sinh^ 




+ iln 


\ cosh 


7r(v^-s"/2) 


1 


\) 


1 cosh 


7r(v^+s"/2) 


sinh^ 




+ C0S2 (f S') 





(25) 



For the symmetric case, 6 = 0, this result agrees with the known transmission probability 
for the Eckart potential 10,3 • 

Fig. 2 shows T as a function of h for fixed a and various energies. Note that for small 
energies the transmission probability can be raised by increasing the asymmetry parameter 
h. While with increasing h the barrier height increases too, the width of the barrier decreases 
leading to a maximum of the transition probability at finite b. 

For e < the solutions of the Schrodinger equation are not bound, except for the case 
when e coincides with an energy eigenvalue. To calculate the energy eigenvalues of the 
potential we seek for poles of the S-matrix at negative energies, or poles of the transmission 
amplitude, respectively. This way we get for the bound state spectrum 

en = -\[-{2n+l) + s']\ (26) 

in which n = 0, 1, . . . , rzmax where nmax is the largest integer with 2n + 1 < s' . From (Q) and 
(pT|) the corresponding bound states are obtained as 

i,^{z) = [t2zY--^' [t2{z - 2^1 {l + n- s\ -n- 1 - |; ^) , (27) 



where z is given by (|T0D and Nn is a normalization constant. Note that the second coefficient 
of the hypergeometric function is a negative integer. Therefore the hypergeometric function 
is a polynomial of z of order n. Of course, ipn is a complex function but one can show that 
the real and the imaginary parts of ipn differ only by a constant. Hence, the eigenvalue is 
not degenerate. Clearly bound states only exist for s' > 1. 
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In the symmetric case 6 = bound states only occur for a < and the spectrum has the 
well known form 

en = -^{-(2n + l) + {l + 4a}^}'. (28) 

For a = the potential (|l|) becomes antisymmetric and the energy eigenvalues of the bound 
states are 



2 



2 



^ -(2n + 1) + + + (46)^1 I . (29) 



.2 2 

In both of these cases there exists at least one bound state for arbitrary a and 6, respectively. 
If a > and 6 7^ the condition for having at least one bound state is 6 > ^/a. 

Fig. 3 shows the transmission probability T as a function of a < for fixed energy and 
various values of the asymmetry parameter h. The transmission probability for 6 = has 
maxima exactly at those points where new bound states are formed, i.e. where e„^^^ is equal 
to zero. Then, for all energies the transmission probability reaches 1 and an incoming wave 
function is completely transmitted. We note that a similar behavior of the transmission 
probability is also found for 6 7^ and a > 0, however, in general the position and the height 
of the maxima depend on the values of e and h. For instance, for e the transmission 
probability has a maximum for a = 1 at 6 = 1 (Fig. 2) where n^ax = 0. Further maxima 
occur for s' = 2?7, + 1, i.e. for a = 1, at those values of b where 

&=^^\/(2^ + 1)^ + 3. (30) 
The transmission probability for s' = 2n + 1 reads in the limit e — > 

2^(0)L'=2n+l = l/cOsh2(^^72). (31) 



B. The Equilibrium Density Matrix 



We assume now that the energy of particles incident from the far left and right is ther- 
mally distributed. If we use the dimensionless formulation introduced in @) and (|^), the 
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state of the system is then described by the equihbrium density matrix at inverse tempera- 
ture (3 which has the coordinate representation 

Pp{x,x') = ^{x\exp{-f3H)\x') (32) 

where Z is an appropriate normahzation constant and H the Hamilton operator from @. 
Using a complete set of orthonormal eigenfunctions of the Hamiltonian with 

de $:(x) ^,{x') + J2 Ki^) ^n{x') = 6{x - x') (33) 



n=0 



one gets 

Pi3{x,x') = ^ 



de exp(-/3e)$:(a;)$,(x')+ E exp(-/3e„)$;(x) $„(a;') 



(34) 

n=0 

Hence, the density matrix at given inverse temperature /3 as a function of x and x' 
is reduced to a sum over the discrete spectrum and an integration about the continuous 
spectrum which both has to be done numerically. 

The bound states in (|34D are given by (^). The scattering states $e are constructed 
appropriately from the solutions \E'(x) in (p!3D. According to ([T9| ) we have for given energy 
two independent scattering states of the system with the asymptotic behavior given by ( |20D 
and (|22D. Hence, t$+/v^27r and t$_/v^27r form an orthonormal set of wave functions for 
e > 0. Thus, according to (P^, the diagonal part of the density matrix, i.e. the position 
distribution P{x) = pp{x,x), reads 



rOO rifr / \ ' •■ma x 

|-exp(-/?A;V2)T(A;)(^|$+(fc,x)|2 + |$_(A;,x)pj + exp(-/3e„ 



TV 

n=0 



where we have introduced the dimensionless momentum k by 



e = \e. (36) 

Performing the integration (^) for various values of a and h one sees that for large values of 
a and 6, corresponding to potential fields with a large height or depth, the numerics becomes 
rather nontrivial because of the strongly oscillating hypergeometric functions. The potential 
(§) has extrema 
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4 V 



fa ± \/a2 + 62 



(37) 



at 



X' 



arsinh 




(38) 



For increasing a and b the absolute values of become large. In dimensional units this 
corresponds to the fact that the extrema of the potential become large compared to the 
quantum mechanical energy scales of the system. The latter are given by the ground state 
energies of harmonic oscillators with frequencies proportional to the curvature of the poten- 
tial in the well and at the barrier top, respectively. Now, for large we have the typical 
situation where a semiclassical approximation should be valid. Thereby, the semiclassical 
expansion in powers of h in dimensional units corresponds to an expansion in powers of 
l/v^2V^ in the dimensionless formulation introduced above. We will show in the following 
that the usual WKB-approximation is indeed sufficient except for a small parameter range 
where the simple semiclassical approximation has to be improved upon. In particular, this 
improved semiclassical result can be handled numerically very simply and provides an excel- 
lent approximation for the exact density matrix in the region where the numerical evaluation 
of Eq. (|35D becomes difficult. 



The dimensionless coordinate representation of the equilibrium density matrix of a quan- 
tum particle moving in the potential V{x) may be written as an imaginary time path integral 



Here, the functional integral is over all paths x(r), < t < (3 with x(0) = x and x{(3) = x'. 
Each path is weighted by its Euclidean action 



III. SEMICLASSICAL DENSITY MATRIX 



A. Path Integral and Semiclassical Approximation 



i 




(39) 
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? n 

S[x] = dr 







+ V{x) 



(40) 



To evaluate the path integral in a semiclassical expansion we first determine the max- 
imum of the weighting factor, that is the minimum of S[x]. This is given by the classical 
action ^[xc/], where Xd is the classical path solving the classical equation of motion following 
from Hamilton's principle 6S[x] = 0. An arbitrary path in ([39| ) reads 



x{t) =x,iiT)+y{T), (41) 

where ?/(r) describes the quantum fluctuations about the classical path. The fluctuations 
have to fulfill periodic boundary conditions y{0) = y{[3) = 0. Using ( PD the full action 
is then expanded around its minimum. This way the dominant part is separated off and 
one is left with a functional integral over periodic paths. In the simple semiclassical or 
WKB-approximation, the expansion of the full action is truncated after the second order 
term leading to an exactly solvable Gaussian path integral 

If there exists a set {x%} of classical trajectories in V{x), this procedure must be per- 
formed for each x";, and all contributions are summed to yield the semiclassical density 
matrix 

Ppix,x')=j:4r^~'^^"'^^ (42) 

V rv 



where = detlS"^ S[x]/ 6x{ti)6x{t2)\x=x'^^} is the determinant describing the Gaussian inte- 
gral over the quantum fluctuations [Q]. Ja is given by the product of eigenvalues {A„} of 
the second order variational operator 6'^ S[x]/ 6x{Tl)6x{T2)\x=x•^^ as 

J„ = 27r/? n(^A°) (43) 

n 

where is an appropriate normalization constant. As long as the second order variational 
operator is positive definite, i.e. A„ > for all n, the Gaussian approximation gives the 
leading order fluctuation term. But a problem arises if one of the eigenvalues A„ tends 
to zero, e.g. as the temperature is lowered. Then, the quantum fluctuations of this mode 
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become arbitrarily large and the simple semiclassical approximation breaks down. Generally, 
the vanishing of an eigenvalue A„ defines a point where new minimal action paths in the 
potential V{x) become possible. This is well-known as the problem of caustics. For our 



purposes we also use another representation of Ja equivalent to (^31). One finds 

J„ = 2vr x^iiP) xS(0) / ——2 dr (44) 

where x denotes the time derivative dx/dr. This way the semiclassical approximation is 
completely determined by the classical paths. 



B. Classical Paths 

In the following we first investigate the classical motion in the inverted potential ~V{x). 
Afterwards, the corresponding classical action and the contribution of the quantum fiuctu- 
ations about the classical path are determined. 

The Euclidean energy of the system reads 

1 9 a + 6sinh(x) , ^ 

The solution of the classical equation of motion for fixed e is given by 



b IP' a 

x{t) = arsinh <( _ + W—^ + _ - 1 sin V2e (r + /5o) (46) 



where Pq is an integration constant. The parameters e and /?o of the solution are determined 
by the boundary conditions x(0) = x, x{[3) = x' leading to 

h rip a 
sinh(x) = ^ + VY^ + ^-1 sm[V2^(3o] 

sinh(x') = ^ + + sin [VTei/S + /3o)] • (47) 

Due to the sine function in (^) one expects that there may be more than one trajectory 
connecting the endpoints for a given inverse temperature or "time" f3. The maximum 
amplitude Xm of these path is then determined by V{xm) = e- Note that e may be negative. 
Therefore /3o is in general an imaginary phase and due to the relation sm{ix) = i sinh(x) the 
oscillatory solution may change into a part of an unbounded motion. 
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1. Symmetric Barrier Potential (b = 0) 



Let us first consider the motion of a classical particle in the inverted symmetric barrier 
potential (a > 0, 6 = 0) |13|. The boundary conditions { ^7\) then reduce to 



sinh(a;) = J — - — - sin y/2e(3^ 



smh{x') = sj'^i^ sin [V2'e{(3 + /5o)] . (48) 

The simplest case is obtained for the boundary condition = where x"*" is the 

coordinate P5| ) at the barrier top. Then, for short times, i.e. high temperatures, there exists 
only the constant solution x(r) = of ( PB| ) with e = = a/2 . When the temperature 



is lowered two new solutions with e = 7r^/2/3^ and /3o = 0,/5 arise for (3 > (3c = 7r/v2y+. 
These solutions describe oscillations in the well of the inverted potential. It will be shown 
below that the constant path, which is stable for high temperatures, becomes unstable for 
times (3 > (3c where the newly emerging trajectories with Xm 7^ are stable (cf. Fig. 4). This 
change of stability is typical for a Hopf bifurcation. 

For X = x' 7^ the bifurcation scenario is perturbed (see Fig. 4). For (3 < (3c one has 
only one stable path that continuously extends also to the region /5 > (ic- For some critical 
value /3c two new branches emerge describing an oscillation to the other side of the well in 
the inverted potential. The solution of smaller amplitude is unstable while the other one is 
stable. The critical value /?c increases with increasing |x| = \x'\. For endpoints within the 
harmonic range of the potential, where |x| <^ 1, the critical temperature /?c remains close to 
(3c. 

Further bifurcations appear in an analog way near /5 = n/3c , n = 1,2,3,.... For x = 
x' = we obtain from ( ^8|) the parameters of the two trajectories bifurcating a.t (3 = n(3c 

b = 0,p/n 

^e(3 = mi n = 1,2,3,.... (49) 

The parameters of the corresponding trajectories for x = x' 7^ describing the perturbed 
bifurcation near /3 = n(ic are determined by 
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sm[V2e{(3 + (3o)] = sm[V2e(3o] 

^/Ye{p + Po) = V2i/3o(-l)" + nil. (50) 

For n odd the newly emerging paths are periodic in coordinate space with x(/3) = a;(0) and 
x{l3) = — x(0), while for n even they are even periodic in phase space with x(/3) = x(0), 
x{P) =x{0). 

In the general case x ^ x' ^ there is again a similar bifurcation sequence as the 
temperature is lowered. However, for high temperatures, i.e. for very small times j3, the 
solution of (48) connecting x and x' is a segment of an unbounded trajectory in the inverted 



potential with e < 0. 

2. Asymmetric Barrier Potential (a,b> 0) 

Now, we investigate the general case a ^ b ^ where we assume a,b > 0. The classical 
motion in the other cases a,b < and a < 0, 6 > as well as a > 0, 6 < is readily obtained 
in an analog way. 

Let us consider the boundary condition x = x' in some detail. Then (^7|) can be trans- 
formed to read 



2e(/? + /5o) = V2e/3o(-l)"+n7r, 

x(0) = X. (51) 

For n odd the phase f3o can easily be determined from ( |5TD and the energy e follows from 



sinh(x) = A + -^2 + Ye i^f 1 (^2) 



For n even the energy e is given by v2e/3 = nn and the phase Pq is determined by 



^^"^ = 2;^ + V w + - ' r ) • 



Now, let us investigate the solutions of ( [5^ ) and ( |5B| ) in detail. For simplicity we first 
discuss the motion near the extrema of the potential [cf. (|37D and (|38|) ]. For x = x~ the only 
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solution of the equation of motion is the constant solution x{t) — x~ with energy e = V~ 
for all temperatures. For x — x'^ we see that for high temperatures there exists only one 
constant solution x{t) = x~^ with e = (see Fig. 5a). When the temperature is lowered two 
new solutions appear at Pc < Pc describing oscillatory solutions to the left side of the well 
in the inverted potential with x{0) < 0. Note that in the symmetric case /3c is always larger 
than Pc- The solution with amplitude Xm closer to x~^ is unstable while the other solution 
is stable. For the stable path \x'^ — Xjn\ increases with decreasing temperature. When we 
arrive the temperature 



the energy e changes sign and the oscillatory solution changes into a segment of an unbounded 
trajectory. For /3 — > oo the amplitude Xm tends to x~ (see Fig. 5a). For the unstable solution 
the difference \x~^ — Xm\ decreases with decreasing temperature until we arrive the critical 
temperature where P — Pc- There the unstable solution and the constant solution meet in 
a double bifurcation point. At P = Pc the constant solution becomes unstable and a new 
solution with growing amplitude appears. This stable path describes an oscillation to the 
right side in the well of the inverted potential. In an analogous way further bifurcations 
emerge with decreasing temperature. 

Now, for X ^ x'^ and x ^ x~ four cases must be distinguished. First, for x~ < x < Xab 
the only paths connecting x and x' with x — x' are segments of unbounded trajectories with 
e < 0. These paths extend from high to low temperatures. 

Second, in the case x > x^ we have at high temperatures only one solution with i;(0) > 
that extends continuously to low temperatures (see Fig. 5a). As the temperature is decreased 
two new solutions with amplitudes x^ < appear at a critical inverse temperature Pc 
describing oscillations to the left side of the well in the inverted potential (see Fig. 5a). One 




(54) 



where Xrr, reaches 



Xab — arsinh(— a/6). 



(55) 



15 



of these solutions is stable and the other one is again unstable. With increasing x the critical 
inverse temperature (3c also increases. At /3 = Pab{x) (Q) the stable oscillation changes to 
a segment of an unbounded trajectory with amplitude Xm that reaches x~ if [3 —>■ oo. The 
amplitude Xm of the unstable path moves away from x+ with increasing (3 until we arrive 
the temperature where | At this point the path that is periodic in 

configuration space [x = x',x{0) = —x{(3)] becomes a path that is periodic in phase space 
[x = x',x{0) = x{(3)] (see Fig. 5a). One can show that this new path and all further paths 
bifurcating with decreasing temperature are unstable. 

Third, we consider the intermediate region Xab < x < x^. (^) defines a minimal inverse 
temperature for which a bifurcation occurs. The corresponding boundary value x is 
denoted by x^. In the symmetric case {b = 0) we find f3^ = (3c and x'^ = x~^ = 0, while in 
the asymmetric case the values < Pc and < x+ depend on the potential parameters 
a and b and have to be calculated numerically. Again two cases must be distinguished. If 
Xab < a; < we have a bifurcation scenario similar to the former case x > x~^. On the 
other hand, in the case x^ < x < x"*", the high temperature solution does not extend to 
low temperatures continuously. At an inverse temperature (3c < (3 < (3c two new solutions 
with Xm < emerge, a stable and an unstable path. When the inverse temperature is 
increased the unstable path and the high temperature solution meet at a bifurcation point 
and disappear while the new stable solution extends to low temperatures. Increasing the 
inverse temperature further another bifurcation of a stable and an unstable solution occurs, 
however, now with x^ > x^. 

Finally, for x < x~ again two cases have to be distinguished (see Fig. 5b). For x near x~ 
the high temperature solutions continuously extend to low temperatures and no bifucations 
occur. As X becomes smaller a bifurcation occurs for x = x~ at a critical inverse temperature 
/9~ > /3c that may determined from (|5^). This scenario remains true for all x < x~. Then 
a stable and an unstable path bifurcate at a critical temperature. As the temperature 
is decreased further the unstable path meets the high temperature solution in a second 
bifurcation point while the stable paths extends to the low temperature region (see Fig. 5b) 
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In the general case x ^ x' the boundary conditions (0) may be evaluated in an analogous 
way and similar bifurcation sequences occur as the temperature is lowered. 



C. Classical Action and Fluctuation Determinant 



Fortunately, not all extremal action paths have to be taken into account in the semiclas- 
sical approximation (E2|) for a given temperature since the corresponding classical actions 



are not always minima of S[x\. One gains from (P0| ) and ( ^61) for the action of the classical 
paths 



Sa = -ef3 + 



2e 



Re 



f 



arctan 



arctan ■ 



/tan V2^(/5 + /5o)/2 +h 



f tan \V2^(3o/^] + h 



where Re denotes the real part. In ([56| ) we have introduced the complex coefficients 



(56) 



/ = i(6 + za) 
b 

9 = 1- + ''' 
4e 



h 



IP' a 



16e2 2e 



- 1, 



(57) 



and the quantities e and /5o are determined by the boundary conditions (^Tf). In the sym- 
metric case (6 = 0) (|56|) reduces to 

' ^ - YePo]]) (58) 



Sci = —eP + y/a ( arctan \ \j — tan v2e(/5 + /3o) 



— arctan { \ — tan 
2e 



Considering first the case a; = a;' = x'^ = we have for (3 < 13c only the trivial solution 
xia) = with the action 



S40]=aP/2. 

At I3c = T^l \fa. two new branches emerge with the action 



(59) 



(60) 
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which for /5 > /3c is smaller than 5'ci[0]. Therefore, as mentioned above, we have a Hopf 
bifurcation at /3 = Pc- The other paths with x = x' = bifurcating at /3 = nPc have the 
action 



n 



(61) 



Hence, the actions of the classical paths appearing at bifurcation points for temperatures 
P > 2/5c differ from S'^j^ by terms which are large compared to 1 in the semiclassical limit 
h 0, which, more precisely, corresponds to the limit —>■ for the present case. As 

will be seen below, also in the asymmetric case at most two classical trajectories contribute 
to the semiclassical approximation of the equilibrium density matrix. 

We now proceed to calculate the Gaussian fluctuations about the extremal action paths. 
After some algebra one obtains from and (EBI) for the determinant J 



J 



2tx 



-1/2 



-1/2 



1 



4e 



+ h sin 



X 



— t: H sm 

8e2 2e / 



b 



1 



4e 



+ h sin 



2e/3o 



/iV2e/3cos V2epo 



V2'ef3] + — /i {cos [v^ePo 
2e(/? + /3o 



— cos 



2e{(3 + (3o)]} 



cos 



(62) 



with h from (|57|) . Using ( ^Tf ) one can show that J is positive for the stable classical paths 
but negative for the unstable one. Of course, for unstable paths the fluctuation determinant 
is only defined as an analytical continuation, or formally by (|4^) and (^31). Further, one sees 
that J vanishes at the inverse temperatures (3 = Pc where the new stable and unstable paths 
emerge. We note that (63) is also valid for e < by virtue of the relations sin(ix) = i sinh(a;) 
and cos{ix) = cosh(x). For instance, for x = x' = x~ , the only solution of (^) is given by 
e = and the fluctuation determinant J reads 



J[x' 



2tt 



2\V- 



sinh J2\V-\P]. 



In particular, when we insert the trivial solution for x 
e = and (|62|) reduces to 



X 



X ' 



(63) 



into ( P7D one gets 



which vanishes for (3 = jSc and is negative for (3 > f3c- 

Let us now discuss the main features of the fluctuation determinant in the symmetric 
case [b = 0). Quahtatively similar reasoning apphes to the general case. For the stable 
branches with x = = bifurcating aX (3 = (3c one finds for P > Pc 

n 

This demonstrates the change of stability near (3 = (3c and shows that the simple semi- 
classical approximation becomes insufficient for temperatures close to (ic = 7r/v2V^, 
i.e. Tc = V2V+h^ /TckBfTiLQ in dimensional units. 

At inverse temperature (3 = 2(3c and for boundary conditions x = again two new 
classical paths bifurcate from the unstable trivial solution x{a) = 0. However, using Eqs. ( ^8[ ) 
and (§2p, one sees that the corresponding determinant J is negative, while J[x'^] becomes 
positive for 27r < y/2V~^(3 < Svr. In function space this means that the paths newly emerging 
aX (3 = 2(3c are unstable in one direction corresponding to one negative eigenvalue of the 
second order variational operator, see (^Sf) , while the trivial solution is unstable in two 
directions with two negative eigenvalues leading to a positive fluctuation determinant. More 
generally, one can show that all classical paths newly emerging for inverse temperatures 
(3 > 2tt/ -\/2F+ are unstable due to the fact that successively the eigenvalues A„ of the 
second order variational operator change from positive to negative values as the temperature 
is lowered. The only stable paths for inverse temperatures (3 > (3c are those bifurcating first 
near /3 = /3c = vr/-\/2V+. These paths are therefore the only ones contributing to the 
functional integral for lower temperatures in the semiclassical limit . 

For temperatures close to Tc a change of stability occurs, as described above. The 
minimal action paths are then not well separated in function space and the semiclassical 
approximation according to (|42|) breaks down. This problem will be addressed in section 4. 
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D. Semiclassical Density Matrix 



Now, for high temperatures, P < and given endpoints x, x' there is only one classical 
path with amplitude Xm which is obtained from (^61) by choosing e and Pq according to (^7|). 
In practice, e and Pq have to be determined numerically. By virtue of Eqs. (|56D and ( |62D 
the semiclassical density matrix for high temperatures, P < P^, is then given by 

Pfs{x, x') = i ^ e-^(--). (66) 

To illustrate this result we have depicted in Fig.^ the diagonal part of the density matrix 
P{x) = pp{x,x) for the symmetric potential with 6 = in the semiclassical approximation 
(|66|) and the exact result ( p5| ) at inverse temperature P = Pc/2 and for two values of a. 
Note that even for a small a = 4, the semiclassical approximation does rather well. The 
probability density P{x) is normalized by the condition lim^^-too P{x) = 1 

For P'^P^ new classical paths may emerge, first for coordinates near the critical coordi- 
nate a;^. Then, one unstable path and two stable ones exist, determined by the solutions e 
and Pq of (^7|) . Apart from a narrow region about the critical inverse temperature where the 



bifurcation occurs, these trajectories are well separated in function space and the sum in (^21) 
contains the contribution of the two stable paths. The density matrix in the semiclassical 
approximation then reads 

1 



Pf3[X,X ) = - 



(67) 



_\Jj{Xml) \Jj{Xm2) 

where we have denoted the amplitudes of the two stable paths with Xmi and Xm2- When 
S{xmi) exceeds S{xm2) by terms that are large compared to 1, it is inconsistent to take both 
contributions in ([67| ) into account. In section 4 we will show that this is the case for all 
coordinates outside a narrow region around the barrier top and around the well region at 
x~ . Hence, mostly the functional integral is dominated by the contribution of one path and 
the semiclassical density matrix reduces to (|66D. 

In Fig.|^ the semiclassical approximation for P{x) calculated from ( p^ ) is depicted to- 
gether with the exact result (^) at inverse temperature P = 2/?c for the potential (3) with 
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a = b = 4. From (pS]) we see that there is only one bound state. Hence, in the semiclassical 
approximation deviations become apparent in the well region, while for the barrier region the 
semiclassical result gives an excellent approximation to the exact result. The semiclassical 
approximation in the well region becomes better if either a or 6 are increased corresponding 
to more bound states in the well. Only for coordinates in the vicinity of the barrier top the 
semiclassical approximation is given by the sum (|67|), which matches for larger x onto the 
simpler formula (|66|). 

As we have seen above, the simple semiclassical approximation breaks down near those 
inverse temperatures (3 where bifurcations of the classical paths occur. As we will see in the 
next section this problem becomes relevant for coordinates in a narrow region around the 
barrier top, near the critical coordinate xjf, and in the well region, near the critical coordinate 
x~ . Only there, the simple semiclassical determination has to be improved due to the fact 
that the relevant classical paths are not well separated in function space. Furthermore, this 



analysis also gives the precise conditions for the validity of the formulae (p6|) and (p7|) 



IV. SEMICLASSICAL DENSITY MATRIX NEAR CAUSTICS 



The specific results presented in Figs. (6) and (7) are for temperatures well above or below 
the critical temperatures (3^ and respectively, where the simple semiclassical results (66) 
and (67) are valid. A more refined treatment is required near the critical temperatures. In 
we have shown how to determine the semiclassical density matrix in the critical region 
near the barrier top. Here, we shall apply these findings to the present problem. 

Near the barrier top the potential (H) with a > may be expanded as 



v{0 = V' 



i-e + la^e + W + o {e) 



3 ^ 2 

where ^ = x — x'^ and the potential parameters are given by 

as = - tanh [x^^ 

4 5 , 2 / 4-\ 

a4 = tanh [x j . 



(68) 



(69) 
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Of course, for coordinates near the barrier top |^| ^ 1, and the anharmonic terms in Eq. 
(|68|) are much smaller than the parabolic one. However, for j3 ^ (3^ they will be found to 
be crucial for the fluctuation integral even in the semiclassical limit. 

To see this more clearly, we use the representation ( ^31) of the determinant J as a product 
of eigenvalues of the second order variational operator. Representing the scaled quantum 
fluctuations in the form 

we have shown in that for high temperatures the second order variational operator is 
transformed to a diagonal matrix with eigenvalues 

C7TTI \ 
Tm-p) '''' 

Then, the path integral over the fluctuation modes reduces to a product of independent 
Gaussian integrals. By virtue of (^) and (|7ID one obtains our earlier result (Q) for J by 
inserting the trivial solution. 

Now, it is obvious that the Gaussian approximation becomes insufficient for inverse 
temperatures near f3c where the eigenvalue Ai vanishes. As the temperature is lowered, 
the fluctuation mode amplitude Yi increases with decreasing Ai. As discussed in [|] the 
anharmonic terms lead to new stability and one gets for inverse temperatures near /3c for 
the action 

S[x] - S[x,i] = 'jJ2 ^nYl + V(em, Yx) (72) 

^ n=2 

with the fluctuation potential 



9 2 

Al — Ac + -04 (^m — C,m,c) 



Here, we have introduced the critical parameters 
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and 

9 

= ^a4^^,c- (75) 
The coefficients 03 and 04 are given by (|69D . In (|73| ) is determined by the cubic equation 

^a4d + l^^sd + Ai^m = ^ + e') • (76) 

This equation has one or three real solutions corresponding to the one or three classical 
paths. In fact, for given endpoints near the barrier top the amplitudes Xm of the classical 
paths determined by (^Tj) and V{xm) = e are given by the roots of Eq. (^) according to 

From Eq. ( [75| ) one sees that the coefficient 

9 2 

Al(Cm) = Ai - Ac + -04 {^rn " ^m,c) (77) 

of the harmonic term may vanish, but the remaining terms of the fluctuation potential 
always constrain Yi to fluctuation amplitudes at most of order (V^Y^^. Introducing the 
sum variable r = (^ + ^')/2 we see from (|76D that for given P the coefficient Ai(^m) vanishes 
at values r = r± (see Fig.|3) where 



'^±(Ai) = Tc 

with 



4 — (-r)" 



(78) 



3a4\/+/3 



2 



C,c- (79) 



Stt 

Note, that Vc = for ^ = On the two curves r±{Xi) bifurcations of classical paths 
with endpoints near the barrier top occur according to the discussion in subsection 3.2.2. 
This analysis gives in detail the regions in the Ai — r— parameter space where one or three 
classical paths exist (see Fig.||). 

For fixed the potential V{C,m,Yi) behaves as follows (see Fig.|^). At high temper- 
atures the fluctuation potential exhibits one minimum at Yi = according to one real 
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solution of ([f6D. When the temperature is lowered to the value where r±{Xi) = r two new 
extrema arise at Yi ^ 0. Then, the fluctuation potential shows two minima and one local 
maximum. In the region in coordinate-temperature space around (rc, Ac) where |Ai(,^m)| is 
smaller than order (\/+)~^/^ both minima differ by terms which are smaller than order 1 
and are not separated by a local maximum exceeding the minima by terms which are larger 
than order 1. 

In this region the functional integral over the fluctuation modes now takes the form p 



V\y\ e(5-5^')[^l 
with the fluctuation integral 



K{U) (80) 



\ 27rsin(V27+/3) 



i^(U = y£ / ^^^^'""^"-^^"^ (81) 



describing the contribution of the marginal mode Y\. The fluctuation integral is flnite for 
positive and negative values of Ai — Ac- 

By virtue of (^6]) and (^O]) the dimensionless density matrix near the top of the barrier 
potential @ and in the vicinity of Ai = Ac takes the form 



KiUe-'^""-^. (82) 



Z ^ 27rsin(v/27+/5) 

Usually, explicit results for the distribution can be calculated only numerically. After solving 
([76|) for given Ai and ^, ^' one has to insert into the fluctuation potential (|73|) and deter- 
mine the fluctuation integral dsTf ). All steps involve only rather simple numerics. Outside 
the critical region where |Ai - Ad < {V+)-^/^ and 1^1, |f | < max( asiV+y^^^ , (V+)-3/8 
) the result ( ^21) smoothly matches onto the various semiclassical approximations specifled 
above in ( |66D and (^). In particular, one can show @] that for coordinates outside the 
critical region, the minimum of the fluctuation potential V{C,m,Yi) corresponding to the 
globally stable path differs from the other minimum by terms which are large compared to 
1. Therefore, for coordinates outside the critical region the density matrix is given by (|66D. 
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For coordinates near the barrier top, ( |66D must be used for Ai — Ac > (y~^) (|67|) for 
Ai - Ac < (V^+)-^/^ and (H) in between. 

In principle, the procedure to determine the fluctuation integral (^) for critical coordi- 
nates in the well region near x~ is similar to the case studied above. As we have seen in 
section 3, a stable and an unstable path newly emerge at the caustic (see Fig. 5b). These 
paths with endpoints near x~ have amplitudes Xm with large \xm — x^\. Hence, the trajecto- 
ries explore the nonlinear range of the well region and the eigenfunctions of the second order 
variational operator cannot be calculated analytically. However, the features of the fluctua- 
tion potential are the same as discussed above for the caustic near x^. Thus the fluctuation 
integral can also be determined in the following way. First, since there is one unstable path, 



the determinant J in (|62D contains only one vanishing eigenvalue at the caustic. Second, 
well beyond the caustic one has one or two stable classical paths which are well separated 
in function space so that the density matrix is given by (|66D and (|67D . Finally, the equation 
for the boundary values (47) expanded around the critical coordinate x~ and the critical 
temperature P~ yields in leading order a cubic equation similar to (81) for = — Xm,c 
where Xm,c is the corresponding amplitude at the critical point. The critical values x~ and 
P~ are determined by (47) and derivatives thereof. Hence, similar to (78), the fluctuation 
potential near the critical point may be written in the form 

V{xm, Y) = h{xm)Y^ + h.Y'' + -^J,Y^ . (83) 

Here, J{xm) denotes the determinant of the second order variational operator of a stable 
classical path with amplitude x^ and Y is proportional to the amplitude of the marginal 
fluctuation mode. The coefficients J3 and J4 are determined by the corresponding cubic 
equation. Then, the density matrix near the caustic is given by 

P/3(a;, x') = exp [-^(a;™)] / ^= exp[-V(a;„, F)], (84) 

JtX) Y ZTT 

where S{xm) is the action (|56|). The regions in temperature-coordinate-space where the 
various semiclassical formulas have to be used may be determined in a way analogous to the 
analysis near the barrier top. 
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Fig. 10 shows the divergent simple semiclassical result (|66D, the improved formula (p2D, 
and the exact result ( p5D for the position distribution P{x) at temperature (3 = j3^. The 
exact result and the improved formula differ by terms which are at most of order 1 / y/2V^. 

The formulas (p2D and (Q) for the improved semiclassical approximation combines with 
the formulas (66) and (67) for the usual semiclassical approximation to yield the semiclassical 
equilibrium density matrix for all temperatures. In particular, the position distribution for 
very low temperatures will be discussed in the next section. 



V. SEMICLASSICAL BOUND STATES AND RESONANCES 

From (34) one sees that for low temperatures the density matrix is dominated by the 
contributions of bound states while the contribution of the extended states is exponentially 
suppressed. For large (3 we expect large deviation between the semiclassical approximation 
(67) and the exact result (35) due to the fact that differences between the exact and semiclas- 
sical energy eigenvalues of the bound states are exponentially amplified in the corresponding 
Boltzmann factors for f3 oo. However, we will show in the following that the semiclas- 
sical position distribution can nevertheless be used to determine the absolute square of the 
wave functions |\I^„(x)p, = 0, l-.n^ax of the bound states with energies En,n = 0, X-.n^ax 
reasonably well. 

For very low temperatures the position distribution function is, apart from exponentially 
small corrections, dominated by the ground state. Thus, the low temperature limit enables 
us to extract the absolute square of the ground state wave function from the shape of 
P^[x). Additionally, from the low temperature behavior of the maximum of Pfi{x) in the 
well region we can extract the semiclassical energy eigenvalue of the ground state -Eq by use 
of the formula 



En = - — lim In 

° A/5 /3-00 



(85) 



where Xmax denotes the coordinate where Pi3{x) has its maximum. The wave function 
|\l/o(x)p is then given by exp(/5o-Eo)-P/3o(^) where Pq has to be choosen sufficiently large so 
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that this product becomes independent of temperature |T^. This procedure can be extended 
also to the next exited state in the following way. We substract the contribution of the 
ground state from the probability distribution at a higher temperature jSi where jSi has to 
be choosen large enough so that the reduced position distribution 

Pi,^(x) = Pp{x) - \^o{x)\^exp{-f3Eo) (86) 

is dominated by the first exited state. From the reduced probability distribution we ex- 
tract the first exited state |\E'i(x)p and its energy Ei. Accordingly, we get a sequence of 
temperatures Po > Pi > (32 > ■■■ and reduced probabilities Pn^pi^x) from which we can find 
the semiclassical approximation for the wave functions |\E'„(x)p and the energies En. To 
illustrate this, we depict in Fig. 11a the exact (27) and semiclassical ground state wave 
function |\E'o(x)p for an antisymmetric potential with potential parameters a = and 6 = 4. 
For these parameters there is only one bound state (cf. (26)) with the semiclassical energy 
Eq = —0.326.., while the exact energy eigenvalue from (29) is e = —0.459.. . The tempera- 
ture choosen to define the ground state is (3q = 12. Note that the wave function is very well 
reproduced, while the Boltzmann factors differ largely: exp(/?o(-E'o — eo)) — 5. 

The semiclassical density matrix can also be used to study scattering states, in particular 
resonances. As we have seen in the previous sections the semiclassical approximation of the 
equilibrium density matrix does rather well even for systems that have only a few bound 
states except for very low temperatures where the amplification of errors in the energy 
eigenvalues has to be taken into account. From the exact solution of the equilibrium density 
matrix we know that the contribution of the scattering states exhibits several peaks. These 
resonances correspond to energy "eigenstates" with complex eigenvalues. In Fig. lib we 
depict the exact and semiclassical position distribution (upper curves) as well as the reduced 
position distribution Pi^p{x) (|86D (lower curves) for the same set of parameters as in Fig. 11a 
at temperature /3 = 0.5. Since for these parameters there is only one bound state, Pi^is^x) 
represents the contribution of the scattering states to the position distribution. According 
to the single bound state in the potential, the distribution Pi^f3{x) shows one resonance 
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peak which is a typical quantum mechanical effect of interference of scattering states. The 
semiclassical approximation, based on purely imaginary time mechanics, reproduces the 
resonance very well. 
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FIGURES 



FIG. 1. The potential (3) for various potential parameters a and b. 

FIG. 2. The transmission probability (6) as a function of b for fixed potential parameter a = 1 
and various energies e. 

FIG. 3. The transmission probability (6) as a function of a at fixed energy e = 0.015 for various 
potential parameters b. 

FIG. 4. Bifurcation diagram for the amplitude Xm of the classical path x(r) versus the 
dimensionless inverse temperature /3 for potential parameters a = 1, 6 = 0. The solid lines are the 
stable and the dashed lines are the unstable branches for a pure {x = x' = 0, thick lines) and a 
perturbed (x = x' > 0, thin lines) bifurcation. 

FIG. 5. The asymmetric bifurcation scenario for various boundary conditions and potential 
parameters a = 1, 6 = 5. a) The thick lines represents the amplitude x^ for x = x^, the thin lines 
represent the amplitude Xm for x > x+. The solid (dashed) lines represent the stable (unstable) 
paths, b) The amplitude x^ for the boundary values (from above) X ^ X ^ x^ , X — x^ and 

FIG. 6. Position distribution P{x) for a symmetric barrier potential (b = 0) at inverse 
temperature with /? = /3c/2 for (a) a = 4 and (b) a = 25. The solid lines represent the exact 
result (^) and the dotted lines the semiclassical formula 



FIG. 7. Position distribution P{x) at inverse temperature (3 = 2(3c- The potential parameters 
are a = 6 = 4. The solid line represents the exact result (|35|), the dashed lines the semiclassical 
formula (|6^). 



FIG. 8. The Ai-r-plane is divided by the two curves r±{\i) (solid lines) into two regions in 



which the cubic equation (|76|) has one or three (shaded region) solutions. Along the line re(Ai) 
one has irn = £.m,c- 
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FIG. 9. The fluctuation potential ( |73| ) for = Cm,c and inverse temperatures near f3c- 

FIG. 10. Position distribution P{x) for a symmetric barrier potential with a = 25 and 6 = 
at the critical inverse temperature = vr/S. The solid line represents the exact result (35), the 
dashed line the divergent semiclassical formula (66), and the dotted line the improved semiclassical 
result (90). 

FIG. 11. a) The absolute square of the ground state wave function for an antisymmetric 
barrier potential with a = and b = 4. b) The position distribution P(x) (upper curves) and the 
resonance (lower curves) for the same set of potential parameters at inverse temperature f3 = 0.5. 
The solid lines represents the exact results, the dashed lines the semiclassical ones. 
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Figure 1, F. J. Weiper et al. 
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Figure 2, F. J. Weiper et al. 
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Figure 4, F. J. Weiper et al. 
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Figure 5a, F. J. Weiper et al. 
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Figure 5b, F. J. Weiper et al. 
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Figure 6a, F. J. Weiper et al. 
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Figure 6b, F. J. Weiper et al. 
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Figure 7, F. J. Weiper et al. 
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Figure 8, F. J. Weiper et al. 
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Figure 9, F. J. Weiper et al. 
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Figure 10, F. J. Weiper et al. 
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Figure 11a, F. J. Weiper et al. 
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Figure lib, F. J. Weiper et al. 
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